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D^- D^ Mixing Analyses at BaBAR 
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INTRODUCTION 

Understanding D meson (charm) mixing is an impor- 
tant step in measuring CP violation in the charm sector. 
It also fills in a gap between the well-measured cases of K 
[1] and B [2, 3] system mixing, both of which have down- 
type quarks in the intermediate state, where charm mix- 
ing has up-type quarks. Since mixing in the D" system is 
expected to be small in the Standard Model [4] (modulo 
the hard-to-predict effects of long-distance interactions 
[5]), charm mixing also offers a chance to observe New 
Physics either through CP violation in mixing [10] or a 
large mass difference between the D mass eigenstates [5] . 
In this proceeding, we summarise the result of four dif- 
ferent approaches to measuring the D mixing parameters 
at BABAR, involving the decays D" -> K+ tt" [6], D" ^ 
K+ K- or 7r+ tt" [7], D^ -^ K+ tt" 7r° [8], and D° -^ 



K+ 



I- tt"!" tt 



[9]. 



DETECTOR 



We present analyses of e+e~ collisions at a center-of- 
mass (CM) energy of 10.58, collected at the BABAR de- 
tector at the PEP-II storage ring. Particle identification 
is done by dE/dx measurements from two tracking de- 
tectors and from measuring Cherenkov angles in a ring- 
imaging detector. D mesons are tagged by reconstructing 
_D*+ "> Z?° 7r+ and D*^ -^ D^ tt^ decays, and assigning 
flavour according to the charge of the slow pion. 



FORMALISM AND NOTATION 

D mesons are produced in pure flavour eigenstates, 
jD") or |Z?°). These flavour eigenstates are not equal to 
the mass and lifetime eigenstates 



\D,) = p\D^> 



q\D") 



\D2) = p\D")-q\D'') 



by which they propagate and decay. Therefore, a particle 
produced as a D'^ may become a £)° before its decay. The 
process is governed by the mass and lifetime differences 
of the Di and D2 states; these decay according to 

|£>i(i)) = e-^("i~*ri/2)*|Di) 
\D2{t)) = e^*("^-*r^/2)*|D2) 

where TOj, Fj are the mass and width of the Di state. We 
define 
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The quantities x and y are collectively referred to as mix- 
ing parameters. Estimates within the Standard Model 
vary from 10^"' (counting only short-distance effects) to 
as high as 1%. Establishing the presence of New Physics 
requires either x » y, or CP violation [5]. 



EXPERIMENTAL APPROACH 

The studies considered here use a common apparatus 
for tagging D mesons as either £)° or D^, and for mea- 
suring their decay times. In particular, by considering 
only D mesons from D* ^ D° tts, we can use the charge 
of the slow pion tts to determine the production flavour 
of the £)", and measure its flight length from the decay 
vertices of the D* and D^ particles. We make use of the 
mass of D^ candidates (todo) and the mass difference 
Am between D° and D* candidates to extract our signal 
yields, and to define sidebands for background studies. 
Figure 1 shows distributions of these quantities for the 
D^ -^ K^ 7r+ analysis, which may be considered typical. 

For historical reasons, D mesons whose decay flavour 
matches their production flavour (e.g. D*^ -^ D^ 7r+ 
with D° -^ K^ 7r+) are called 'right-sign' (RS), while 
the opposite case is referred to as 'wrong-sign' (WS). 



Wrong-sign decays may come about either through mix- 
ing or through doubly-Cabibbo-suppressed (DCS) Feyn- 
mann diagrams. To distinguish the two cases, we use the 
decay-time distribution, as wiU be shown for each decay 
mode. 

In addition to these two sources of wrong-sign events, 
there is the case where a correctly reconstructed D° is 
matched with a pion not from a D* decay to produce 
a spurious D*\ this is referred to as the "mistag" back- 
ground. Another source of background is D mesons re- 
constructed with the correct tracks, but wrong particle 
assigments, or with tracks missing; this is the "bad D°" 
or "mis-reconstructed charm" background. Finally there 
is background from combinatorics. 

D° -^ K+ TT- 

In the limit of small mixing and CP conservation, the 
decay-time distribution for wrong-sign decays of mesons 
produced as D^ may be approximated as 

^^^ oc Ru+y'VR^{Tt) + -{x'^+y'^)(Ttf{l) 
e ^ ^ 4 

where x' and y' are related to x and y by 

x' = X cos 5ktc + y sin 5k-k 
y' ~ y cos Sj^^ — x sin S^tt ■ 

The angle Sk-w is the strong phase between Cabibbo- 
favoured (CF) and DCS decays. The quantity Rjj is 
the amplitude, in the absence of mixing, for the D^ to 
decay by a DCS process; the term quadratic in t is the 
amplitude, in the absence of DCS processes, for the D" 
to mix and then decay as a -D"; and the term linear in t 
is the interference term between these two processes. 

We apply Equation 1 in two ways: The first is to en- 
force CP conservation by fitting both D° and D'^ sam- 
ples together. The second is to search for CP violation 
by doing two fits, calculating a;'^ and y' for D^ and Z?" 
separately. 

We use 384 fb~^ of e+e~ data, pairing tracks of op- 
posite charge to make Z?" candidates, and then pairing 
these with slow pion tracks to make D* candidates. The 
phase space available for slow pions is small; we require 
their momentum to be greater than 0.1 GcV/c in the 
lab frame, and less than 0.45 GeV/c in the CM frame. 
We fit the full decay chain, constraining the D* to come 
from the beam spot, the D'^ and slow pion to come from 
a common vertex, and the K^ and tt to come from a 
different common vertex. We reject candidates if the x^ 
probability of this fit is less than 0.1%. The D° decay 
time and error on the decay time are taken from this fit; 
candidates whose decay-time error exceeds 0.5 ps are as- 
sumed to be badly reconstructed, and thrown away, and 
we also require that the decay time be between —2 and 



4 ps. We further require the CM momentum of I?° can- 
didates to be at least 2.5 GeV/c, which suppresses back- 
grounds from B-meson decays and combinatorics. Where 
multiple D* candidates share tracks, we use only the can- 
didate with the highest x^ probability from the fit. With 
these criteria, our samples consist of 1,229,000 RS and 
64,000 WS D° and D° candidates. Figure 1 shows their 
distribution in rriKTr and Am. 
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FIG. 1: a) tuktv for wrong-sign (WS) candidates with 
0.1445 < Am < 0.1465 GeV/c^ and b) Am for WS candi- 
dates with 1.843 < TTiK^ < 1.883 GeV/c^ The fitted PDFs 
are overlaid. 



We extract the mixing parameters using an unbinned, 
extended maximum- likelihood fit, which proceeds in 
three stages. The first step is to fit the rriKTr — Am dis- 
tributions to extract shape parameters in these variables; 
these are then fixed in subsequent fits. Next we fit the 
RS sample to extract the D^ lifetime and resolution func- 
tions, using the ttix^ — Ar?i parameters from the previous 
step to separate the components. Finally we fit the WS 
sample for the mixing parameters using three different 
models. The first model assumes no CP violation and no 
mixing; the second permits mixing, but not CP viola- 
tion; the third allows both mixing and CP violation. 

The rriKTr — Am distributions are fitted to a sum of 
four PDFs, one each for signal, mistags, bad 2?° and 
combinatorial background. Of these, the signal peaks in 
both m^K-K and Am. The mistagged events - correctly 
reconstructed D° with a pion not from a. D* decay - 
peak in m^T^ but not in Am. Bad D'^ events have a 
D'^ with one or more daughters missing, or assigned the 
wrong particle hypothesis; they peak in Am, but not in 
niK-n- Finally, combinatorial background does not peak 
in either variable. Figure 1 shows these various shapes. 
The signal peak contains 1, 141, 500 ± 1, 200 candidates 
for the RS sample, and 4, 030 ± 90 for the WS. 

We describe the decay-time distribution of the RS sig- 
nal with an exponential convolved with a sum of three 
Gaussians, whose widths are proportional to the mea- 
sured event-by-event error on the decay time. The com- 
binatorial background is described by a sum of two Gaus- 
sians, one of which has a power-law tail; the mistag back- 
ground is described by the same PDF as the signal, be- 
cause the slow pion has little influence on the vertex fit. 



For the WS signal, we use Equation 1, convolved with the 
resolution function determined by the RS fit. Figure 2 
shows the data, overlaid by these various PDFs. From in- 
spection, it is clear that the fit allowing mixing describes 
the data better than the one which imposes zero mixing. 
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FIG. 2: a) Projections of the proper-time distribution of com- 
bined D" and D" WS candidates and fit result integrated over 
the signal region 1.843 < m^vr < 1.883 GeV/c^ and 0.1445 
< Am < 0.1465 GeV/c^ The result of the fit allowing (not 
allowing) mixing but not CP violation is overlaid as a solid 
(dashed) curve, b) The points represent the difference be- 
tween the data and the no-mixing fit. The solid curve shows 
the difference between fits with and without mixing. The dif- 
ference between the mixing- allowed fit and the data is there- 
fore the difference between the solid curve and the points, or 
essentially zero. 




FIG. 3: The central value (point) and confidence-level (CL) 
contours for 1 - CL = 0.317 (la), 4.55 x 10"^ {2a), 2.70 x 
10"^ (3cr), 6.33 X 10"^ (4cr) and 5.73 x 10"'^ {5a), calculated 
from the change in the value of —2\nC compared with its 
value at the minimum. Systematic uncertainties are included. 
The no- mixing point is shown as a plus sign (-I-). 



TABLE I: Results from the different fits. The first uncertainty 
listed is statistical and the second systematic. 



Fit type 


Parameter 


Fit Results (/lO 


-') 


No CP viol, or mixing Rd 


3.53 ± 


0.08 ± 


0.04 


No CP 
violation 


Rd 


3.03 ± 


0.16 ± 


0.10 


/2 

X 


-0.22 ± 


0.30 ± 


0.21 


f 

y 


9.7 ± 


4.4 ± 


3.1 




Rd 


3.03 ± 


0.16 ± 


0.10 


CP 


Ad 


-21 ± 


52 ± 


15 


violation 


X'2+ 


-0.24 ± 


0.43 ± 


0.30 


allowed 


y'+ 


9.8 ± 


6.4 ± 


4.5 




x'2- 


-0.20 ± 


0.41 ± 


0.29 




y'- 


9.6 ± 


6.1 ± 


4.3 



Figure 3 shows the likelihood contours of the mix- 
ing parameters from the fit allowing mixing but not 
CP violation, including systematic uncertainties. The 
point of maximum probability is in the unphysical re- 
gion where x'^ is negative; adjusting for this by moving 
to the most likely point in the physical region, .t'^ = 0, 
y' ^ 6.4 x 10~3, we find that -2Aln/: is 23.2 units be- 
tween the most likely physical point and the point of 
no mixing. Including the systematic uncertainties, we 
thus find mixing at a significance of 3.9 <t. Table I 
shows the results of our fits in more detail; we find no 
evidence for CP violation, as shown by the asymmetry 
Ad = [Rd - R'd)I{Rd + ^d) where subscript '+' indi- 
cates only the D^ sample was used, and '— ' indicates the 
D'^ sample. 

We evaluate systematic uncertainties froin three 
sources: Variations in the fit model, in the selection cri- 
teria, and in our procedure for dealing with track-sharing 
D* candidates. The most significant source of systematic 
uncertainty in Rjj and the mixing parameters is from 
the fit model for the long-lived background component 



caused by other D decays in the signal region, followed 
by the presence of a non-zero mean in the time-resolution 
function, caused by small misalignments in the detector. 
For the asymmetry Ao, the dominant contribution is un- 
certainty in modeling the differences between K'^ and 
K~ absorption in the detector. 



D° -^ K+K- OR 77+ 77" 

For D mesons decaying to CP cigenstates, mixing 
changes the decay time distribution in such a way that 
we may, to a good approximation, consider the decays 
exponential with changed lifetimes ([11]) 

T+ = T^^ [I + \q/p\{y con (j)f — xsa\(j)})\~ 
T^ = T° [I + \p / q\ {y con 4) f + X sin 4) f)]~ 

where t° is the lifetime for decays to final states which 
are not CP cigenstates, and r"*" (r~) is the lifetime for 



D° (£)") decays to CP-even states. Wc can combine the 
three Ufetimes into quantities 

ycp = t°/{t) - 1 

AF = (r°A.)/(r). 

Here 0/ is the CP-violating phase 0/ = aig{qAf /pAf)^ 
Af {Af) being the amphtude for D" (-0°) decaying to the 
final state /. (r) is the average of t+ and t~ , and Ar is 
their asymmetry (t+ — t^) / (t+ + t^). In the absence 
of mixing, both ycp and AY are zero. In the absence of 
CP violation in the interference of mixing and decay (ie, 
(pf = 0), Ay is zero and ycp = U- 

For this analysis, we use 384 fb~^ of BaBAR data, and 
measure the lifetimes for the CP-even decays^ _D" — > A'+ 
K~ and I?" ^ 7r+ tt^, and for _D° -^ A'" 7r+, which is 
not a CP eigenstate and thus gives our r". 

In addition to particle identification requirements, the 
cosine of the helicity angle (defined as the angle between 
the momentum of the positively charged D*^ daughter in 
the Z?° rest frame, and the D^'s momentum in the lab 
frame) is required less than 0.7; this suppresses combina- 
torial backgrounds. D° candidates are then combined 
with pions to produce D* candidates. Electrons are 
rejected by combining pion candidates with each other 
track in the event and vetoing those which form a good 
photon conversion or pion Dalitz decay, as well as by 
dE/dx measurements. The requirements for slow pions 
and the vertex fit of the D* are the same as for the 
£)" -^ K~ 7r+ analysis (Section 5). 
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FIG. 4: Reconstructed D^ mass distributions for the three 
D" samples, within ±0.8 MeV/c^ of the At7i peak. 

Figure 4 shows the mass distributions of D^ candi- 
dates; Table II shows the yield and purity of the samples, 
calculated using events within a 15 MeV/c^ D^ mass and 



Sample Size Purity (%) 
7r+ 730,880 
-K+ 69,696 
■K--K+ 30,679 



K- 
K- 



99.9 
99.6 
98.0 



^ Charge conjugation 
noted. 



is implied throughout unless otherwise 



TABLE II: Sample sizes and purities. 



0.8 MeV/c^ Am window. We fit the decay time distri- 
butions of these samples using an unbinned maximum 
likelihood fit to all five decay modes simultaneously, us- 
ing separate PDFs for signal decays, mistagged events, 
mis-reconstructed charm events, and combinatorial back- 
ground. 

As with the D^ -^ K^ 7r+ study, we model the decay- 
time distribution of signal events using a simple exponen- 
tial convolved with a sum of three Gaussians for the res- 
olution. Each Gaussian has a width proportional to the 
event-by-event error on the measured decay time; their 
mean is common, and allowed to be offset from zero to 
account for any effects of detector mis-alignment. Mis- 
tagged events - that is, events with a correctly recon- 
structed I?*^, but wrongly assigned slow pion - account 
for about 0.4% of the sample; of these, half will have 
the wrong flavour assignment to the D^ . However, they 
have the same decay-time distribution and resolution as 
true signal. Hence we model these events using the signal 
PDF, but reversing the flavour assignment. 

Mis-reconstructed charm events have an exponential 
decay-time distribution, which we convolve with a sin- 
gle Gaussian. The fraction of such events is obtained 
from simulation, which we check by comparing data and 
Monte Carlo in the sidebands 1.89 < mj^o < 1.92GeV/c^ 
and 0.151 < Am < 0.159 GeV/c^. We estimate the charm 
background as (0.009 ± 0.002)% of events in the signal 
region for D° -^ R- 7r+, (0.2 ± 0.1)% for D° -^ K+ 
K-, and (0.15 ± 0.15)% for D° -^ 7r+ tt". For combi- 
natorial background, we model the decay-time distribu- 
tion as the sum of a Gaussian and a modified Gaussian 
with a power-law tail, the latter accounting for a long- 
lived component. Each decay mode has its own shape 
for combinatorial background, the shapes being deter- 
mined from flts to the sideband regions; the fraction of 
this background is again estimated from Monte Carlo 
with uncertainties derived from comparison of MC and 
data. We find (0.032 ± 0.003)% in the D° -^ R- 7r+ 
mode, (0.16±0.02)%in D^ -^ K+ A'", and (1.8±0.2)% 
inD° -^ 7r+ tt". 

We consider several sources of systematic error, in- 
cluding variations of the signal and background models, 
changes to the event selection, and detector effects. We 
vary the models by changing the signal PDF shape and 
size, as well as the position of the signal box. We also 
test our resolution model by forcing the common mean 
of the three Gaussians to zero, and by allowing it to float 
separately for different bins of the D^ polar angle. Of 



TABLE III: Summary of systematic uncertainties on ycp and 
AY, separately for KK and tt+tt^ and averaged over the two 
CP modes, in percent. 





o 


yc. (%) 


CTAY (%) 


Systematic 


KK 


TT+TT Avg. 


KK 


TT+TT Avg. 


Signal model 


0.130 


0.059 0.085 


0.072 


0.265 0.062 


Charm bkg 


0.062 


0.037 0.043 


0.001 


0.002 0.001 


Comb, bkg 


0.019 


0.142 0.045 


0.001 


0.005 0.002 


Selection criteria 


0.068 


0.178 0.046 


0.083 


0.172 0.011 


Detector model 


0.064 


0.080 0.064 


0.054 


0.040 0.054 


Quadrature sum 


0.172 


0.251 0.132 


0.122 


0.318 0.083 



these effects, the largest systematic uncertainty derives 
from widening the D'^ mass window, which increases the 
amount of badly-reconstructed signal events in the sam- 
ple. 

We vary the mis-reconstructed charm model by chang- 
ing its fraction in the fit, by varying its effective lifetime, 
by using a different sideband region, and by using a decay 
time distribution obtained from Monte Carlo instead of 
the sideband data. Due to the purity of the data, these 
effects are all small, the largest being from varying the 
background fraction in the D*^ -^ 7r+ tt^ mode, where 
the purity is worst. 

We vary our event selection criteria in two ways: By 
throwing out or keeping all multiple candidates (as op- 
posed to selecting the candidate with the best x^ proba- 
bility for its vertex fit), and by changing the acceptable 
range of errors on decay times. The last, which changes 
the amount of poorly reconstructed signal events, has the 
largest effect. 

Finally, we consider effects of our understanding of 
the detector by repeating our analysis with different mis- 
alignment parameters. This changes our fitted lifetimes 
by up to 3 fs; but since the lifetimes change by simi- 
lar amounts, and we are considering ratios of lifetimes, 
the effect on the mixing parameters is small. All these 
systematic effects are summarised in Table III. 

The results of these decay-time fits are shown in Ta- 
ble IV. From the measured lifetimes, we extract 

ycp = 1.24±0.39(stat)±0.13(syst)]% 
AY = [-0.26±0.36(stat)±0.08(syst)]% 



Mode 


Lifetime (fs) 


D" ^ K- n+ 


409.33 ± 0.70 


JJO p*+-) _^ ^+ ^- 


401.28 ± 2.47 


D° (LI*-) ^ K+ K- 


404.47 ± 2.52 


JJO (J3,*+) _^ ^+ ^- 


407.64 ± 3.68 


L>° (D*-) ->7r+ -K- 


407.26 ± 3.73 



which is evidence for D^-D^ mixing at the 3-sigma level, 
and consistent with CP conservation. This amount of 
D^-D^ mixing is consistent with Standard Model predic- 
tions. 



Z?" 



K^ 



For the case of D^ decays to three-body final states, we 
can modify Equation 1 to give a decay-time distribution 
for each point in the decay phase space: 



A{P,t) 



-,-rt 



\Ar 



(2) 



TABLE IV; Measured lifetimes for the different decay modes. 
Uncertainties are statistical only. 



+ \Ap\\Ap\{y" cos 5p ~ x" sin 5p ) Tt 



+ \Ap\'{x"^+y"^){Tty 



In analogy with Equation 1, Ap is the amplitude (in the 
absence of mixing) for D^ mesons to decay by a DCS 
process to the point P on the Dalitz plot. The term 
quadratic in time is the amplitude (in the absence of 
DCS processes) for the D° to mix before its decay, and 
then decay to the point D by a CF process. Within this 
term, the factor Ap is the amplitude for the CF decay, 
while the remaining factors are the mixing amplitude. 
The term linear in time is the interference between the 
DCS and mixing terms. The quantity 5p is the phase 
of the intermediate states in the decay, relative to some 
reference resonance. As with the D^ — > K~ 7r+ case, 
an unknown strong phase Sxtt-^o between CF and DCS 
decays prevents us measuring x and y directly; instead 
we are sensitive to 

x" = X cos 5k-ktt° + y sin Skutt" 
y" = y cos 5ktt-k» - x sin Sk^-^" ■ 

As with the previous two analyses, we use 384 fb~^ 
of BaBAR data, reconstructing D^ —^ K~ tt"*" tt" candi- 
dates from two oppositely-charged tracks and two photon 
candidates with energy at least 100 MeV. The 7r° can- 
didate is required to have a lab momentum of at least 
350 MeV/c, and a mass-constrained fit probability of at 
least 1%. The slow pion is required to have a momentum 
transverse to the beam axis of at least 120 MeV/c, and 
the D^ candidate to have a CM momentum of at least 2.4 
GeV/c. As in the previous two analyses, we extract the 
-D" decay time, with error, from a vertex fit constraining 
the D* to the beam spot; this fit is required to have a x^ 
probability of at least 1%. 

Figure 5 shows the rriKn-n and Am distributions that 
result from these criteria. We fit these distributions as 
described for the D^ -^ K^ 7r+ study in Section 5; the 
fit to the WS sample uses shape parameters from the 
RS fit, suppressing the associated systematics. Table V 
shows the yields for each component. 

We compute the quantity Ap in Equation 2, the time- 
independent amplitude of CF decays to the point P on 
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FIG. 5: (Distributions of RS (top) and WS (bottom) data 
(points with error bars) with fitted PDFs (dashed line) over- 
laid. The mxTTTT distribution (left) requires 0.145 < IS.m < 
0.146 GeV/c^; the Am distribution (right) requires 1.85 < 
rriKTVTv < 1.88gevcc. The white regions represent signal events, 
the light gray (blue) misassociated tt^ events, the medium 
gray (red) correctly associated tt^ with misreconstructed D" 
events, and the dark gray (green) remaining combinatorial 
background. 



Category 


N events (RS) 


N events (WS) 


Signal 


639802 ± 1538 


1483 ± 56 


Combinatoric 


1537 ± 57 


499 ± 57 


Mistag 


2384 ± 57 


765 ± 29 


Bad D" 


3117 ±93 


227 ± 75 



TABLE V: Number of RS and WS events of signal and back- 
ground in the m^o and Am signal region. 



the Dalitz plot, by fitting the RS Dalitz plot to an iso- 
bar model, using the signal and background fractions 
obtained in the fit to the rriK-m^ ~ Ajti distribution. 
The background PDF is empirically determined from the 
TOjf TTTT ~ A?Ti sidebands, and its fraction is set to the back- 
ground fraction derived from the rriKTi-n — Am fit. 

With Ap (or more accurately, the phases and ampli- 
tudes for intermediate resonances from which Ap can be 
calculated) known, we then go on to fit the WS sample 
simultaneously to the Dalitz plot. We thereby determine 



FIG. 6: WS D" decay-time distribution (crosses) with fit 
(solid blue line) overlaid. The green and red regions show 
the mistag and (combinatoric-t-bad-D") backgrounds respec- 
tively. These backgrounds are taken from sideband data, 
which accounts for their jagged shape. 

Ap, and the decay-time distribution, to extract the mix- 
ing parameters. The signal decay-time PDF is taken as 
Equation 2 convolved with a sum of three Gaussians, 
as described for the previous two analyses; the parame- 
ters of the Gaussians are extracted from a fit to the RS 
decay-time distribution, and fixed in the WS fit. For the 
background components, mistagged events are described 
by the RS parameters, since they contain correctly re- 
constructed D° mesons; the other two background com- 
ponents are described empirically using the sidebands. 
Figure 6 shows the WS fit projected to the decay time. 
From this we extract x" = 2.39 ± 0.61 (stat.) ± 0.32 
(syst.)% and y" = -0.14 ± 0.60 (stat.) ± 0.40 (syst.) 
%. This excludes the no-mixing hypothesis at the 99% 
confidence level. 



D° 



K+ 



As in the case of D" -^ K^ 7r+ tt*^, the four-body final 
state K^ -K^ 7r+ tt^ has a decay-time distribution which 
varies across the phase space. However, since the phase 
space is enlarged by one dimension, and the data sample 
for this analysis is smaller, we do not fit for a decay time 
at each phase-space point. Instead we integrate across 
phase space to get a WS to RS decay-rate ratio of 



- WS 



(t) 



rRs(t) 



Rd + ay' ^ RoiTt) 



+ t{i'^ + y'^){Ti)^ 



(3) 



where a tilde indicates integration over phase space. The 
quantity a is a suppression factor accounting for strong- 
phase variation across phase space; in effect it measures 
the amount of information we lose by the integration pro- 
cedure. As in the D^ -^ K~ tt"*" analysis, Rd is the am- 
plitude for doubly-Cabibbo-suppressed decays, the term 



quadratic in time is the amplitude for mixed decays, and 
the term hnear in time is the interfcrence between the 
two. Again we account for an unknown strong phase by 
using variables 

i! = X cos (5 + 2/ sin 5 
y' = y cos 5 ^ X sin 5 

where 5 is the strong phase difference integrated across 
phase space. Equation 3 assumes CP conservation. To 
account for possible CP violation in interference between 
DCS and mixed contributions, we introduce the inte- 
grated CP-violation phase (^, and parametrise CP vio- 
lation in the mixing itself with \p/q\. This allows us to 
make the substitutions 



ay 

,,2 



b/g|±i (ay' 



cos 4) ± fix' sin t 



{x^+y^) - \p/qf^{x'+y^) 

in Equation 3, applying plus signs for the D^ sample and 
minus signs for D'^ . (3 is an information-loss parameter 
analogous to a, in this case accounting for phase-space 
variation in (p. 

This analysis uses a 230.4 fb^^ BaBml dataset. The 
reconstruction procedure is analogous to that of the pre- 
vious three analyses, the main difference being the re- 
quirement that neither pion pair have an invariant mass 
within 20 MeV/c^ of the K^ mass of 0.4977 GeV/c^. We 
demand a I?° CM momentum requirement of at least 2.4 
GeV/c. Two vertex fits are performed, one for the D^ 
candidate, required to have a x^ probability of at least 
0.5%, and one for the full D* decay tree. For the latter, 
from which we derive our decay-time value and error, the 
D* is constrained to come from the beam-spot, and the 
probability is required to be at least 1%. The mean at 
for signal events is 0.29 ps; events with Uf > 0.5 ps are 
rejected. The signal yields are calculated from a fit to the 
(77i/f37r, Am) distribution; Table VI shows the results. 



TABLE VI: Signal yields determined by the two-dimensional 
fit to the (rriKSTv , Am) distributions for the WS and RS sam- 
ples. Uncertainties are calculated from the fit. 
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FIG. 7: Distributions of WS data with fitted PDF overlaid. 
The light gray shows mis-reconstructed charm, the medium 
gray shows mistagged events, and the dark gray shows com- 
binatorial background. 
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FIG. 8: Left: Contours of constant Aln£ = 1.15, 3.0 in terms 
of the doubly Cabibbo-suppressed amplitude and the time- 
integrated mixing rate. Right: Contours of constant A In i3 = 
1.15, 3.0 in terms of the normalized interference term and the 
integrated mixing rate, for the D° and D" samples separately. 
The hatched regions are physically forbidden. 



D" 



D" 



WS (1.162 ±0.053) X 10'' (1.040 ± 0.051) x 10^ 
RS (3.511 ± 0.006) X 10^ (3.492 ± 0.006) x 10^^ 



The (mi<-3^, Atti) fit which extracts the signal yields 
also determines shape parameters for those two variables; 
these are then used in a three-dimensional fit which also 
includes the time distribution. The decay time func- 
tion for RS events is a simple exponential convolved 
with a double Gaussian, with widths proportional to 
at and separate means. For mistagged events we use 
the RS decay-time PDF; for mis-reconstructed D° com- 
ponent we use the signal PDF; and for combinatorial 



background a Gaussian with a power-law tail. We fit 
the RS sample to determine the D'^ lifetime and the 
time-resolution parameters, which are then held fixed in 
the fit to the WS sample. We allow yields and back- 
ground shape parameters to vary. Figure 7 shows the 
WS decay-time distribution and fit. Figure 8 shows con- 
tours of constant likelihood in the {Rd^ Rm) plane; we 
find Rm - (0.019+°°^^ 



_Q Q]^5 ± 0.002)% assuming CP conser- 



vation, and Rm = (0.017to;oi6 ± 0.003)% with CP vio- 
lation allowed. There is no significant difference between 
the D^ and Z?" samples in the CP-allowed fit. 

To extract a consistency with the no- mixing hypothesis 
from these contours is not quite straightforward, because 



Mode 



Luminosity [fb 



Mixing 



CP violation 



D^ ^K- 


7V + 


384 


D° ^ K- 


K+ or 7r+ TT" 


384 


D° -* K- 


^+^« 


384 


D" ^ K- 


^+ T,- t:+ 


230.4 



3.9 a No evidence 

3.0 a No evidence 

Exclude NM at 99% CL No evidence 

Consistent with NM at 4.3% CL No evidence 



TABLE VIL Summary of results. 



the linear term in Equation 3 becomes unconstrained as 
Rm approaches zero. We therefore estimate the con- 
sistency of our data with no-mixing using a frcquentist 
method; we generate 1000 data sets of 76300 events each, 
setting the mixing parameters to zero in the generation. 
We then apply our fit procedure to these sets; in 43 cases 
we find an Rm equal to or greater than for the data. 
We therefore conclude that our data are consistent with 
no-mixing only at the 4.3% confidence level. 

We investigate systematic uncertainties from four 
sources, listed in order of decreasing significance. First 
is the at threshold, which we increase from 0.5 to 0.6 ps. 
Second is the decay-time resolution function; we change 
this by fixing one of the Gaussian widths to be exactly 
equal to at , letting the other constant of proportionality 
float as before. Third, the rriKs-n distribution of the back- 
ground is changed from exponential to a second-order 
polynomial. And fourth, we use the nominal value of the 
D^ lifetime instead of the one obtained from our RS fit. 
Taken all together, these uncertainties are smaller than 
the statistical uncertainty by a factor of five. 



SUMMARY AND OUTLOOK 

BABAR has found evidence for mixing in several chan- 
nels, as summarised in Table VIL With the total BABAR 
luminosity expected to reach 750 fb~^ before shutdown, 
or nearly twice the largest amount used in these studies, 
we expect to be able to improve these measurements of 
the D mixing parameters, and to add other channels as 
well. 
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